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ABSTRACT 

We calculate and discuss the light element freeze-out nucleosynthesis in high 
entropy winds and fireballs for broad ranges of entropy-per-baryon, dynamic 
timescales characterizing relativistic expansion, and neutron-to-proton ratios. 
With conditions characteristic of Gamma Ray Bursts (GRBs) we find that deu- 
terium production can be prodigious, with final abundance values 2 H/H > 2%, 
depending on the fireball isospin, late time dynamics, and the effects of neutron 
decoupling-induced high energy non-thermal nuclear reactions. This implies that 
there potentially could be detectable local enhancements in the deuterium abun- 
dance associated with GRB events. 

Subject headings: gamma rays:bursts — nucleosynthesis — stars:neutron 

1. Introduction 

In this paper we show that the high entropy relativistic flows thought to accompany 
catastrophic stellar endpoint events like binary neutron star mergers (Eichler et al. 1989) 
and collapsars (MacFadyen & Woosley 1999) may be a source of interesting light element 
synthesis. In particular, we find that under the right circumstances a significant fraction 
of the relativistic ejecta can be turned into 2 H. Though this cannot be a cosmologically 
significant 2 H source on account of the small mass and presumed low event rates associated 
with, e.g., GRBs, it could represent an appreciable local enhancement. Because the collision- 
less shocking (and subsequent synchrotron emission) of ultra-relativistic winds is thought to 
give rise to observed GRBs, deuterium production in these winds raises the possibility of a 
nucleosynthetic signature of the GRB environment. 
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The general nucleosynthesis associated with rapid adiabatic expansion and freeze-out 
from Nuclear Statistical Equilibrium (NSE) at high entropy-per-baryon, s, was first consid- 
ered in the landmark paper by Wagoner, Fowler, & Hoyle (1967) (hereafter WFH). Those 
authors concentrated on the environments associated with exploding supermassive objects, 
where s/kj, ~ 10 3 , and on Big Bang Nucleosynthesis (BBN), where s/ki, « 10 10 . In both 
supermassive objects and BBN the conditions are expected to be proton-rich (preponderance 
of protons over neutrons) and the characteristic expansion timescale large (e.g., r dyn ~ 100s 
for BBN). The physics of relativistic outflows and potential nucleosynthesis in these sites has 
also been the subject of some recent attention (Pruet, Fuller & Cardall 2001; Thompson, 
Burrows, & Meyer 2001; Otsuki et al. 2000). 

Here we extend the WFH study. We consider freeze-out from NSE over a wide range 
of entropy-per-baryon spanning that in WFH, and a range of neutron-to-proton ratios, all 
the way from proton-rich to neutron-rich. NSE freezeout in these scenarios is calculated 
for expansion timescales ranging from those appropriate for relativistic flows from compact 
objects to those associated with BBN (10~ 6 s < r dyn < 100s). 

We find that the general nuclear physics features of BBN are a recurring theme through- 
out all of these parameter ranges. In particular, 2 H synthesis can be significant for relativistic 
flows. The expected possibly neutron-rich conditions of a GRB fireball are, in a sense, the 
isospin mirror of those for BBN. However, we also show that it is not sufficient to simply 
consider only the NSE freeze-out nucleosynthesis in some parameter regimes appropriate 
for, e.g., GRB fireballs. In fact, the dynamics of these fireballs can differ dramatically on a 
microscopic scale from the conditions treated by WFH. 

In particular, for initially neutron-rich material, the few protons accelerate with the 
relativistic photon and e ± pair fluid to a very high Lorentz factor (7 ~ 10 2 — 10 3 ). In this 
scenario, there will be high energy (~ GeV) collisions of protons on "left-behind" neutron 
stragglers (Fuller, Pruet, & Abazajian 2000). It has been shown that high energy collisions of 
this kind can result in significant destruction and/or production of light nuclei (Dimopoulos 
et al. 1988). We show below that these non-thermal nuclear reactions can modify signifi- 
cantly the simple BBN-like NSE freeze-out abundances of 2 H and other species. In the right 
conditions, 2 H number fraction yields can approach Yd ~ 10%. This is a staggeringly high 
yield given the known, and small, primordial 2 H abundance. 

The deuteron is a particularly interesting nucleus because it has a binding energy of only 
p»2.2 MeV and is notoriously fragile and difficult to synthesize. At present it is thought that 
essentially all of the deuterium in the universe is primordial. Estimates of the primordial 
2 H/H come from measurements of the 82 km/s isotope shift in the Lyman a line in a handful 
of high redshift Lyman limit clouds. These are clouds that intersect our lines of sight to 
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distant QSO's. The values of 2 H/H inferred from these clouds are consistent with each other 
and average to 2 H/H ~ 3 • 10~ 5 (O'Meara et al. 2001). Because the deuterium yield in the 
big bang is a sensitive function of the entropy per baryon (or alternatively baryon to photon 
ratio 77), the observed 2 H/H places important constraints on conditions during the epoch of 
BBN, and gives a determination of the baryon rest mass contribution to closure (scaled by 
the Hubble parameter in units of 100 kms _1 Mpc _1 ), Vt^h 2 pa 0.02. 

Estimates of the present day interstellar medium values of 2 H/H in different directions 
from the UV Lyman a transition are also consistent with each other and average to about 
2 H/H pa 1.5 ■ 10~ 5 (McCullough 1992), in accord with the idea that stellar processing results 
in a net destruction of deuterium. There have also been a number of careful radio searches 
for the hyperfine transition in deuterium. So far these have only yielded upper bounds 
about a factor of four higher than the abundance measured by UV satellites. Nonethless, 
radio searches for pockets of 2 H in supernova remnants may represent a promising avenue for 
detecting evidence of relativistic flows. Ultimately whether or not deuterium from relativistic 
flows is detectable depends on the interaction of the flow with the ISM. We do not calculate 
this here. 



2. Nucleosynthesis in relativistic flows 

The evolution of a hot plasma is governed by the initial radius i?o of the plasma, the 
initial temperature of the plasma T , and the entropy per baryon s/kb in units of Boltz- 
mann's constant. Initially the plasma accelerates under the influence of thermal pressure. 
For optically thick flows the acceleration ceases when the thermal energy of the fireball is 
converted to kinetic energy. Energy conservation gives the final Lorentz factor 7^ of the 
fireball as 

7/ = E /M. (1) 

Here E is the total initial energy (thermal plus rest mass) in the fireball and M is the mass 
of the baryons confined to the fireball (in this and in what follows we set H — c — 1). In 
this paper we will focus on the case where the initial plasma is baryon poor, so 7/ > 10. 
Non-relativistic or mildly relativistic outflows with modest entropy per baryon (s/kb ~ 100) 
can give rise to interesting heavy element synthesis (Meyer et al. 1992; Woosley et al. 
1994; Takahashi, Witti, & Janka 1994), but are relatively uninteresting as far as A<8 nuclei 
are concerned. In terms of 7/ and the initial temperature, the entropy per baryon in the 
fireball is 
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Here m p is the proton rest mass. At these entropies, and for To > 0.2MeV, nuclear 
statistical equilibrium determines that initially the fireball consists of free neutrons and pro- 
tons. The neutron to proton ratio in these relativistic ejecta remains an open and important 
question because it may be a probe of conditions near the fireball source and so help pin down 
the central engine. When weak processes are rapid compared to the expansion timescale, 
n/p will come to equilibrium with these processes at n/p~ 1 for T> a few MeV. The rate of 
electron(positron) capture on free protons (neutrons) is ~ 5(T/lMeV) 5 s _1 . This means that 
e ± capture will influence n/p if (T/lMeV) 5 > 10 4 (10 6 cm/i2 ). (Below we will identify R 
with the dynamic timescale.) This result implies that weak processes may or may not influ- 
ence n/p depending on the details of the fireball source. For relativistic winds from NS-NS 
mergers, for example, neutrino (and e^) processes are too slow to influence n/p, so that the 
composition of the outflow mirrors the composition of the material ablated/ejected from the 
neutron stars (Pruet, Fuller & Cardall 2001). Rough estimates of n/p from different stellar 
and compact object sources range from n/p 10 for NS-NS mergers to n/p ss for flows 
from supermassive object collapse. At any rate, this is one of the few environments where 
free neutrons may be present. We will treat n/p as an unknown free parameter in this work. 

The acceleration of the initially hot plasma has been investigated by a number of authors. 
To a good approximation the Lorentz factor and temperature in the plasma evolve simply as 
7 = R/Rq = Tq/T during the acceleration stage of the fireball's evolution and as 7 = 7/ for 
R > Rolf (Paczynski 1990). These scaling relations imply that an observer comoving with 
the fireball sees the temperature evolve as T = T exp(—tR ) . This determines the dynamic 
timescale as Td yn = Rq 1 - Compact objects typically have radii of order 10 6 cm implying 
Tdyn ~ 10~ 5 s. The dynamic timescale for flows from jets breaking out of massive stars 
during stellar collapse may be several orders of magnitude larger (MacFadyen & Woosley 
1999). The jet evolution in these collapsars is too complicated to allow reliable statements 
about nucleosynthesis. The validity of a simple picture for the jet evolution and associated 
nucleosynthesis depends on the extent of the interaction between the outgoing jet and the 
stellar envelope. 

As the relativistic flow expands and cools, NSE begins to favor bound nuclei. As long 
as the reactions responsible for keeping nuclei in NSE are fast compared to the expansion 
timescale, the nuclear abundances will track statistical equilibrium abundances. Once the 
reaction rates become slow compared to Td yn , the material freezes out of NSE and breaks into 
successively smaller clusters of nuclei in Quasi Statistical Equilibrium (QSE). Below we will 
calculate freeze-out nucleosynthesis yields in relativistic fireballs using the code developed 
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by Kawano for BBN (Kawano 1992). First, however, we illustrate some of the salient aspects 
of nucleosynthesis in these relativistic fireballs by using a simplified toy reaction network. 

Suppose for a moment that deuterium were the only possible bound nuclear species 
(i.e. flow to heavier nuclei is prevented). In this case the evolution of the deuterium number 
fraction Yd (Yx is defined as the number density of nucleus X divided by the number density 
of baryons) is given by 



= PbY p Y n {av) pQ -Dj - Y^Y D \ y « ((Tv) pnDl (pbYpY n - 5 • 10 9 Tg /2 exp(-Q/T 9 )Y D ) (3) 

Here (a"i>) pn D 7 ~ 4.7 • 10 4 cm 3 /(smol) is the reaction rate for n(p, 7)d, T 9 is the tem- 
perature in units of 10 9 K, and Q ~ 2.2 MeV is the binding energy of the deuteron. To 
a first approximation, Yd tracks the equilibrium value, Y^ SE , until the destruction rate 
A 7 Y 7 ~ 5 ■ I0 9 exp(—Q/Tg)(av) pil £> 1 falls below the expansion rate. (The expression for 
the reverse rate follows from the formalism developed by Fowler, Caughlan, & Zimmerman 
(1967).) At this time, Y D freezes out at approximately Y^° ps Y p Y n Tg(T^ 5 /s 5 ). Here r_ 5 is 
the dynamic timescale in units of 10~ 5 s and S5 is the entropy per baryon in units of 10 5 . This 
result implies that for (r_ 5 /s 5 ) > 1CT 1 , i.e. for essentially all reasonable fireball parameters, 
there is time to form interesting abundances of 2 H. In this case the final 2 H abundance is 
limited by the flow to heavier nuclei. 

To estimate the effects of flow to heavier nuclei we employ the code developed by Wag- 
oner (WFH) and subsequently modified by Kawano. This code contains a reaction net- 
work for 26 light nuclei and rates for the relevant weak processes. The parameters of the 
"Universe" relevant for BBN are: Hubble constant H=r dy * l , a freely chosen initial n/p (in 
contrast to BBN where n/p is initially determined by weak reactions and assumed lepton 
numbers/spectra), and a final baryon to photon ratio 77 ~ 3.6/(s/kb)- Note that the 77 rele- 
vant for these relativistic winds is some five orders of magnitude larger than the primordial 
value (77 ~ 6 • 1CT 10 ), and that the dynamic timescale for relativistic winds can be seven 
orders of magnitude smaller than the dynamic timescale of ~ 100s characterizing the epoch 
of BBN. 

In Figure 1 we show the evolution with time of the mass fractions of 4 light nuclei 
for winds with an entropy per baryon s/k b = 10 5 and differing dynamic timescales r dyn = 
10~ 3 , 10~ 4 , 10" 5 and 10~ 6 s. This figure illustrates the effects of the competition between 
production and photo-dissociation of 2 H for fast outflows, and the flow to heavier nuclei for 
slow outflows. For Td yn = 10~ 6 s, the particle capture reactions out of deuterium are too 
slow to influence 2 H/H, and Y D freezes out at approximately 1/80. For r dyn = 10~ 5 s, 
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flow to heavier nuclei, principally via d(d, n) 3 He, d(d, p)t, and t(d, n) 4 He (here t represents 
the triton, d represents the deuteron, and the three rates make roughly equal contributions) 
becomes important and depletes about 90% of the deuterium produced through p(n, 7)d. Of 
these depleted deuterons, approximately 0.1% are lost via d(p, 7) 3 He (or d(n, 7)t for neutron 
rich material). The evolution is similar for the slower outflows, with flow to heavier nuclei 
becoming increasingly important with time. The reaction d(p,7) 3 He begins to dominate as 
the final Y D decreases and the ratio Y p Y D a[d(p, 7) 3 He]/YDYDcr[d(d, n) 3 He] « 10~ 5 Y P /Y D (at 
T 9 = 1) approaches unity. This is the case for the Td yn = 10~ 3 s wind and near this Td yn the 
final Yd is very sensitive to small changes in the entropy and dynamic timescale. This is also 
the regime most like BBN. 

In Figure 2 we display the final 2 H abundance for the range of entropy and dynamic 
timescale of relevance to astrophysical sites. For this figure we assume an initial electron 
fraction Y e = 0.5. Note that for expansion timescales comparable to the free neutron lifetime 
the final yields are sensitive to the temperature at which the neutron to proton ratio is 
specified. We show these very long expansion timescales only for illustrative purposes. The 
two upper contours in Figure 2 correspond to deuterium production being limited by d(7, n)p, 
and the two lower contours correspond to deuterium production being limited by flow out 
of 2 H to heavier nuclei principally by the strong interactions. 

The thick dashed line in the upper half of the graph delimits the region where dynamic 
neutron decoupling occurs. Dynamic neutron decoupling occurs because photons and elec- 
tron/positron pairs dominate the dynamics of the accelerating flow. Protons and charged 
nuclei are coupled to the relativistic light particles electromagnetically, but free neutrons 
are uncharged and have a small magnetic moment, so they are dragged along only by n-p 
collisions. When these collisions become slower than the dynamic expansion timescale, the 
neutrons are left to coast while the protons accelerate without them. In the absence of bound 
nuclei, the condition that dynamic neutron decoupling occurs is that T dyn n p (crv ) np < 1 before 
the end of the acceleration phase of the fireball's evolution (here n p is the number density of 
protons). Putting in the scaling relations describing the evolution of the fireball gives this 
condition as s/10 5 > 3(r dyn F p /10~ 5 ) 1 / 4 (e.g., Fuller, Pruet, & Abazajian (2000)). 

The precise neutron decoupling condition depends on the nuclear composition of the 
flow. If, for example, all of the nuclei are locked into a particles, then neutron decoupling 
cannot occur. Roughly, the fraction of nucleons bound in nuclei is small (< 10%) for Td yn < 
10~ 4 s and entropies along the neutron decoupling line, and the fraction of nucleons bound 
into a particles approaches saturation for r dyn > 10~ 4 and entropies along the neutron 
decoupling line. This implies that for longer dynamic timescales neutron decoupling will 
only occur if the flow is neutron rich, or if the entropy is a factor of a few higher than that 
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Fig. 1. — Evolution with time of the light element mass fractions for four different wind 
cases. For each case the entropy per baryon is 10 5 and Y e is 0.2. 
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given by the neutron decoupling line in the figure. 

If the flow is neutron rich, neutron decoupling will result in high energy (several hundred 
MeV) neutron-nucleus collisions which will destroy and synthesize nuclei. These non-thermal 
reactions are not described by the Kawano reaction network. The net result of these collisions 
depends on the nuclear composition at the time of decoupling. 

We will consider two instructive limiting cases: (i) deuteron rich and alpha poor, and 
(ii) deuteron poor and alpha rich. Case (i) occurs for short dynamic timescales and for 
entropies well above the neutron decoupling line, while case (ii) occurs for longer dynamic 
timescales and entropies near the neutron decoupling line. The final nucleosynthesis is also 
sensitive to the (unknown) details of neutron decoupling. In particular, we will see that 
the final deuterium abundance is sensitive to the exact number of high energy interactions 
suffered by the average nucleus. This number is certainly between one and a few, but the 
exact number can be obtained only with a detailed transport calculation. 

Any deuterons present when neutron decoupling occurs will be destroyed. Likewise, any 
a particles present will be broken apart. Because the branching ratio for a+(~ GeV nucleon) — > 
2 H + X is high, about 50% (Dimopoulos et al. 1988), the production of 2 H through the spal- 
lation of a's may (over) compensate the loss of deuterium in direct collisions. Certainly for 
case (i), the a poor case, neutron decoupling will result in a net decrease in deuterium. For 
case (ii), where a particles dominate, neutron decoupling will drive Y D to approximately half 
of the initial (i.e. at the beginning of neutron decoupling) value of Y4 Ue , provided that the 
average nucleus only undergoes one high energy collision. If the average nucleus undergoes a 
few high energy collisions, the nucleosynthesis is somewhere inbetween the one collision case 
and the fixed point case discussed by Dimopoulos et al. (1988). 

For example, if initially the material is composed entirely of a particles and free neutrons, 
and each a particle suffers one high energy collision with a neutron, then ~ 10% of the 
material will be converted to 2 H. Regions in Figure 2 above the decoupling line where the 
freeze-out mass fraction of a particles is large could approximate this case. If additionally 
this initial spallation takes place right at the neutron decoupling point, then none of the 
deuterium produced from the fragmenting of a will be destroyed. Only if these conditions 
are met can the final 2 H yield be as high as ~ 10%. 

Note that here there is no source of significant post-nucleosynthesis photo-dissociation 
of deuterium. This is because on average there will only be ~ lOOMeV per baryon in 
the form of electromagnetic radiation coming from the decay of pions produced in inelastic 
nucleon-nucleon collisions. This is three orders of magnitude smaller than the ~ lOOGeV 
per baryon coming from the decay of massive particles discussed by Dimopoulos et al.. In 
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addition, here the energy injection by tt and [i decay is well approximated as occurring at 
the decoupling temperature, T D « 0.005(s 5 /r_ 5 r e ) 1 / 3 MeV . This means that, unlike the 
case where a massive particle decays with a lifetime of many system dynamic timescales, 
the maximum energy of the nonthermal breakout photons does not keep rising with time. 
Interestingly, the recognition that the branching ratio for 2 H production in high energy 
nucleon-a collisions is so large suggests that relativistic flows might lead to direct deuterium 
production via collisions with dense winds or ejecta surrounding the site of the relativistic 
wind (Fuller, Guiles, k Pruet 2002). 

In Figure 3 we illustrate the effect that changing the neutron to proton ratio has on 
the nucleosynthesis. Two illustrative examples are shown. The upper curve corresponds to 
a high entropy, fast expansion and there is little synthesis of nuclei heavier than 2 H. In this 
case the final deuterium abundance follows a clean, simple freezeout from NSE, i.e., the Saha 
equation: F D oc Y p Y n ks Y e (l — Y e ), where Y e is the number of electrons per baryon. The 
lower line in Figure 3 corresponds to a lower entropy and may be thought of as a convolution 
of the upper curve with the reactions d(d, n) 3 He and d(d, p)t. The effect of destruction is less 
pronounced on the wings of the curve simply because the deuterium abundance is smaller 
and the reactions d(d, n) 3 He (and other comparable reactions) are slower. 

At this point we can address in broad stroke the issue of Deuterium production in 
proposed GRB models. For long bursts, observed redshifts imply a total energy output in 
gamma rays in GRBs of order 10 53 — 10 54 (fi/47r)erg. Here Q is the opening angle subtended 
by the relativistic outflow. There is some evidence that long GRBs are typically beamed 
by a factor of Q/4tt 1/500 and that the energy output in gamma rays in long GRBs is 
burst-independent and about 10 51 erg. For short GRBs no redshifts have been measured and 
little is known about the energy scale. Not all of the energy E can be converted to gamma 
rays and E is larger than the energy in gamma rays by at least a factor of 5-10. An estimate 
of the energy scale gives the initial temperature in GRBs 



The Lorentz factors invoked to explain GRBs lie in the range ~ 10 2 — 10 3 , with the 
smaller values corresponding to internal shock models for long bursts and the larger values 
corresponding to external shock models for short bursts. This implies entropies per baryon 
in the range 10 4 — 10 6 . From Figure 2 we see that GRBs, and relativistic flows in general, 
associated with compact objects will be accompanied by substantial deuterium production. 
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Fig. 2. — Final deuterium yield as a function of the dynamic timescale and entropy per 
baryon characterizing the relativistic flow for Y e = 0.5. The shaded contours were calculated 
under the assumption that the neutrons remain well coupled kinematically to the flow. The 
upper two lines correspond to 2 H production being limited by photo-dissociation while the 
lower two lines correspond to 2 H production being limited by the flow to 4 He. The thick 
dashed line is the neutron decoupling line as described in the text. Above this line the 
assumption of kinematically coupled neutrons breaks down and high energy nucleon-nucleus 
collisions will change the nucleosynthesis. 
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Fig. 3. — Final deuterium yield as a function of the electron fraction. The upper curve is 
for r dyn = lCT 5 s and s/k b = 10 5 , while the lower curve is for r dyn = lCT 5 s and s/k b = 10 4 . 
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3. Conclusion and Discussion 

There are three length scales discussed in connection with GRBs. The first is the size 
of the central engine (Rq ~ 10 6 cm for compact objects, Rq ~ 10 11 cm for collapsars). This 
region has been the focus of many detailed attempts at understanding the dynamics of the 
collapsing object giving rise to the relativistic flow. The second length scale is the scale 
associated with the shocking that generates the observed signal, 7/-R0 for internal shocks 
or ~ 10 15 cm for external shocks. For obvious reasons, this scale also has been the subject 
of many studies. In this paper we have considered the third, intermediate scale (~ 7-Ro) 
associated with the acceleration of the thermal fireball. This region is rich in nuclear and 
particle phenomena and is the only region sensitive to, e.g., the neutron to proton ratio in 
the flow. 

We have explored nucleosynthesis for the broad range of dynamic timescales, neutron 
to proton ratios, and entropy per baryon which may be found in these relativistic flows 
surrounding stellar endpoint events. For relativistic flows for which the neutrons remain 
kinematically well coupled, the nucleosynthesis is similar to the freeze-out nucleosynthesis 
discussed in connection with BBN. The final mass fraction of deuterium can be greater than 
a few percent, some three orders of magnitude larger than the primordial value. For flows 
for which the assumption of kinematically well coupled neutrons breaks down, light element 
synthesis is dominated by high energy non-thermal nucleon-nucleus collisions. In these latter 
cases the deuterium number fraction can be high, Yd < 10%. This local overabundance of 
deuterium will reside in the shell propagating into the ISM at late times. Recently, these 
shells have gained some attention because a radio identification of a young (~ 100 yr) non- 
spherical HI shell might tell us something about GRBs and their association with supernovae 
(Paczynski 2001). 

Whether or not the deuterium produced in stellar endpoint events is detectable depends 
ultimately on the degree of mixing with the ISM. The total amount of mass in deuterium 
produced in these events is M D « 10~ 6 M Q (£/10 52 erg)(100/7/)(F d /10~ 2 ), where E is the 
total energy in the relativistic flow. For HI column densities of order 10 19 / cm 2 , this deuterium 
will lead to a factor of two increase in 2 H/H over the primordial 2 H/H as long as it mixes 
in a volume less than about a cubic light year. Structures of this size associated with 
nearby supernova remnants are seen and resolved with the VLA. However, even with the 
most optimistic assumption that every supernova is accompanied by a relativistic wind, 
the mechanism discussed here can only contribute to the overall present day deuterium 
abundance at the level of about 0.1 percent. 

It is perhaps frustrating then that such a prodigious yield of deuterium has little leverage 
on the scale of the cosmological abundance. However, the potentially large 2 H production 
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in relativistic flows could possibly be detectable and, therefore, provide a new probe of, e.g., 
GRB central engine physics and expansion dynamics. 
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